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ABSTRACT 

The gamma-rays emitted from the moon or any similar body carry informa- 
tion on the chemical composition of the surface layer. The elements most easily 
measured are K, U, Th and major elements such as O, Si, Mg and Fe. The ex- 
pected fluxes of gamma-ray lines are calculated for four lunar compositions and 
one chondritic chemistry from a consideration of the important emission mec- 
hanisms: natural radioactivity, inelastic scatter, neutron capture, and induced 
radioactivity. The models used for cosmic-ray interactions are those of Reedy 
and Arnold (1972) and Lingenfelter et al (1972), The areal resolution of the 
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experiment is calculated to be around 70-140 km under the conditions of the 
Apollo 15-16 experiments. Finally, a method is described for recovering the 
chemical information from the observed scintillation spectra obtained in these 
experiments. 
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INTRODUCTION 

A planetary body whose atmospheric thickness is not too large emits gamma- 
ray lines whose intensity can be observed by orbital remote sensing and used to 
determine the mean concentration of a number of elements in the upper tens of 
centimeters of the surface being observed. Experiments to gain knowledge about 
the chemical composition of the Moon, Mars, Mercuiy, and other bodies in this 
way have been proposed many times. Up to the present only the Moon has actu- 
ally been studied. Metzger et al (1964) placed an experiment on Rangers 3, 4, 
and 5. None of these was able to make measurements near the moon, but the 
deep space flux was successfully recorded. Vinogradov et al (1967) placed a 
small detector on Luna 10, which orbited the moon. Successful experiments 
were carried out on Apollo 15 and. 16 (Metzger et al , 1973). 

In the present paper we calculate the expected spectrum of gamma-ray line 
emission from the moon, as observed from a detector above the lunar surface, 
for a range of known and possible chemical compositions, primarily with ref- 
erence to the Apollo missions. 

The experimental technique, and the present theoretical treatment, are 
equally applicable to any atmosphereless planetary body such as Mercury, an 
asteroid, or a planetary satellite. For the special case of Mars, whose atmos- 
pheric depth is of the order of the photon interaction length, significant correc- 
tions for absorption are required, but the experiment is still quite practical 
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(Metzger and Arnold, 1970; Surkov et al, 1971). The modifications in the theory 
required for planetary bodies other than the moon are not great, and will be 
noted below where appropriate, 

Gorenstein and Gursky (1970) have made a calculation of gamma-ray 
line emission from the moon, and have reviewed the previous liter- 
ature. Armstrong (1972) has used a Monte Carlo method to derive the flux of 
the stronger lines and of the continuum from the moon. The present work is 
based on a more detailed theoretical model developed mainly for lunar sample 
studies by Reedy and Arnold (1972). The reader interested in a more complete 
discussion of this and other models of cosmic -ray interaction with the lunar sur- 
face should consult that paper. 

In this paper we calculate the fluxes for a number of gamma-ray lines pro- 
duced by the-major sources of line radiation in the moon: decay of naturally 
occurring primordial radioactive nuclei (K, U, and Th) , prompt processes (in- 
elastic scattering and neutron capture) involving secondary neutrons produced 
by galactic-cosmic-ray bombardment, and decay of radionuclides produced by 
galactic-cosmic-ray and solar-cosmic-ray bombardments. The fluxes from a 
semi-infinite plane are calculated for five representative chemical compositions: 
Apollo 11, Apollo 12, Apollo 14, lunar anorthosite, and chondrite. The spatial 
resolutions and fluxes relative to an isotropic detector in orbit about the moon 
are determined for the major sources of gamma-rays. 
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We emphasize the line spectrum because it contains the chemical informa- 
tion. However, it is now clear that the continuum emitted from the lunar sur- 
face has a substantially larger flux than the sum of all lines in the energy region 
of interest (Metzger et al, 1973). Thus it must be accurately determined in 
order to derive correct line intensities. 

Once the production mechanisms are understood, and the relation between 
line intensity and chemical abundance defined as well as possible by theory, there 
remains the problem of deriving the actual photon flux from the observed pulse 
height spectrum in the detector. We give a systematic account of the methods 
used to obtain these fluxes (including the determination of the continuum, treated 
here as an unknown) and hence derive chemical composition, from the experi- 
mental data. 

Sources of Line Radiation. The most direct observations of chemical com- 
position which are possible by this technique are those of the elements K, Th, 
and U. These elements and their radioactive daughters emit gamma-rays as 
given in Table 1. These "natural’' gamma-rays are emitted at a steady and pre- 
dictable rate. Except for disequilibrium caused by radon emanation which at the 
present writing seems to be relatively unimportant for gamma emission processes 
(Gorenstein and Bjorkholm, 1972), the calculation of emission rates is quite 
simple. 
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The other important sources of line radiation all originate in the bombard- 
ment by high energy particles (cosmic -rays). Because of the high energy of the 
cosmic-rays, every possible mechanism of nuclear excitation is active. In 
particular, the positron annihilation line at 0.511 MeV is excited by a wide 
variety of mechanisms. For this reason little chemical information can be de- 
rived from its intensity. In other cases only a small number of elements can 
produce a given line . 

High energy particles from outside the solar system, the so-called galactic- 
cosmic-rays (GCR), bombard the lunar surface continuously, with intensity 
varying slowly over the solar cycle. The interactions of this primary flux with 
a solid body give rise to numerous secondary particles. Those of greatest im- 
portance to gamma-ray production are the particles in the energy range 1-100 
MeV, which consist almost entirely of neutrons. The spectra of the primary and 
secondary pai'ticles are discussed in detail by Reedy and Arnold (1972). 

Of the many possible mechanisms for emission of nuclear gamma-ray lines, 
only three are important for our purposes. These are neutron inelastic scat- 
tering, neutron capture, and induced radioactivity. Neutrons with energies above 
the first excited state of a nucleus can inelastically scatter, exciting states which 
promptly de-excite by gamma-ray emission. There are also several cases of 
prompt gamma- rays from other neutron-induced (n,X 7 ) reactions. Low energy 
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neutrons can be captured by nuclei and produce gamma-rays via (n, 7 ) reactions. 
Nuclear reactions can also produce radioactive species, some of which emit 
gamma-ray lines in high yields. 

The sun also emits energetic particles, the so-called solar-cosmic -rays 
(SCR). These are also chiefly protons, usually having energies of less than 
100 MeV. The flux is intermittent, a given storm lasting only a few days. 

These storms are concentrated in the years around solar maximum. Only the 
gamma-rays resulting from the decay of induced-radioactive nuclei are calcu- 
lated here. 

Chemistry . The gamma-ray experiment is capable of measuring the com- 
position of a lunar or planetary surface over a very wide range. To be realistic 
we have chosen in this paper to present calculations for five chemical composi- 
tions. The analyses performed at the Surveyor V site (Turkevich et al, 1969) 
are so similar to the more detailed data on the Apollo 11 soil that it is reasonable 
to accept the latter as representative of Mare Tranquillitatis, over which the 
Apollo 15 Command Service Module flew. Thus the Apollo 11 results represent probable 
"ground truth". The Apollo 12 soil gives the best information we have on the composi- 
tion of the area of Oceanus Procellarum. The chemical compositions used for 
the Apollo 11 and 12 lunar soils are the average values compiled by J. Warner, 

MSC, as given by Schnetzler and Nava (1971). The Apollo 14 soil (Schnetzer and 
Nava, 1971) from the Fra Mauro region, corresponds closely to what Hubbard 
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and Gast (1971) have called KREEP. They suggest that this type of material may 
be widespread on the lunar surface. The "lunar anorthosite" component has 
been suggested by Wood et al. (1970) as a major constituent of the highlands. 

Finally, we have used an ordinary chondritic composition (Mason, 1962) 
very different from anything yet encountered or expected as a main component 
on the moon. We consider it as a possible asteroidal or planetary surface 
composition. The chemical compositions are given in Table 2. 

Calculation of the Total Emitted Flux . Before discussing the expected 
gamma-ray source strengths as a function of depth, we will give the method of 
calculating the flux emitted for a given source strength. We are not here con- 
cerned with scattered photons, which contribute only to the continuum, but only 
those which undergo no interactions before reaching the surface. There are 
two different quantities that can be used to represent the intensity of photons 
leaving the lunar surface— the flux, the intensity of photons reaching a unit- 
area isotropic detector above the lunar surface, and the current or albedo, the 
intensity of photons passing through a unit-area of (or parallel to) the lunar sur- 
face. Only the flux will be discussed and calculated below. 

In the expressions below, the depth below the surface, x will be in units of 
g/cm 2 and the gamma-ray source strength at depth x is A(x) in units of dis- 
integrations per unit mass and time . The calculations ■will be for 
a monoenergetic photon which passes through lunar material with an exponential 
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mass attenuation coefficient p (in units of cm 2 /g) . All calculations will be for 
a semi-infinite plane source. 

The flux of photons at an isotropic detector above a semi-infinite plane re- 
sulting from photons emitted at a depth x below the surface is 

I(x) = f 2n d<}> f ? sin 6 dd sec 9 e^ x sec 9 , (1) 

J o Jo 4ir 

where A(x)/47 r is the rate of photon emission per unit solid angle at depth x, 9 
is the angle jaoth respect to the normal, and e' MX sec 0 is the fraction of the 
photons emitted from a depth x and with angle 6 that escape from the surface. 
The sec 6 term is the area of the moon's surface that a unit area of detector 
projects for an angle 6 with respect to the normal. Equation (1) integrates over 
solid angles to give the contribution of the flux at the detector from a source at 
depth x in the moon 

I(x) = e 0 Gux). (2) 

e 0 (px) is a Gold integral (Rossi, 1952) and is the same as -Ei(-jux), where Ei 
is the exponential integral. For a source whose activity is constant with depth. 
Equation (2) can be integrated over depth to yield the total flux at the detector, 

1 = v <» 

The expression for the current from a depth x is the same as Equation (1), but 
without the sec 6 term and, when integrated over all angles and depths for a 
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constant source, gives as the current A/4 p, which is half the corresponding 
value for the flux (Equation (3)), 

In the calculations for all sources but natural radioactivity, Equation (2) 
was used to derive the flux at each depth and the integration over depth was 
done numerically. A series (about 50) of planar layers were used, with many 
layers near the surface and fewer layers at greater depths. A(x) is calculated 
for each source for the depth x in the middle of each layer and I(x) determined 
for that layer. I(x) is then multiplied by the thickness of the layer and these 
products for all layers are summed to get the total flux from the semi -infinite 
plane. Using a constant source in this numerical integration results in a 
value for the flux within one percent of the exact value from Equation (3). 

The values used for the mass attenuation coefficient p were those for alu- 
minum (White, 1952), which is a very close parallel to that for natural silicate 
rocks. In the region of interest from 0. 5 to 10 MeV p decreases from about 
0.083 cm 2 /g (mean thickness of 12 g/cm 2 ) to 0.023 cm 2 /g (mean thickness of 
about 44 g/cm 2 ). For a source uniform with depth, half the flux at the detector 
comes from depths of less than 0.27 of the mean thickness (1/p) and 90% from 
depths of less than 1.3 mean thicknesses. 

CALCULATED FLUXES 

To a first approximation the expected flux of a line due to a particular ele- 
ment might be expected to be proportional to its concentration. There are 
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effects, however, which can cause a departure from a simple proportionality 
to abundance. Some of these are small in realistic cases. Thus the range of 
density and mean atomic number is small enough in probable examples that the 
development of the nuclear cascade, and the mans absorption coefficients of the 
emitted gamma-rays, are not significantly affected. However, an iron-nickel 
object, or one made up of "ices," would require a new analysis. The most im- 
portant non-linear phenomena threated here involve the thermalization and 
capture of neutrons and are discussed below. 

Natural Radioactivity. The naturally occurring primordial radioactive 
nuclei _ ura nium, thorium, and their daughters, and K 40 — were assumed to be 
uni formly disturbed to a depth of several meters at the lunar surface. All 
daughter isotopes were assumed to be in equilibrium with their parent, and the 
radon daughters of uranium and thorium were assumed not to have diffused from 
the location at which they were formed. For the emission of photons from a 
uniform source, the flux at an isotropic detector above a semi-infinite plane was 
calculated using Equation (3). The value used for the source strength A in 
Equation {3) for a specific gamma-ray was the disintegration rate of the parent 
nucleus multiplied by the yield of the gamma-ray per disintegration. The im- 
portant parameters and the calculated line intensities for the gamma-rays re- 
sulting from the decay of naturally occurring radioactive nuclei are given in 
Table 1. The calculated fluxes agree quite well with those of Armstrong (1972), 
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and, after dividing by the factor of 2, with the currents calculated by Gorenstein 
and Gursky (1970). The gamma-ray intensities for the decay of V s0 , Lu 176 , and 
La 1 3 8 were found to be negligible. 

Solar-cosmic-ray Induced Radioactivity . The solar-cosmic- rays are emitted 
irregularly from the sun in a few intense solar flares around the period of solar 
maximum during the 11 -year solar cycle. Only those gamma-rays that result 
from the decay of radionuclides produced by solar protons need be considered, 
since high intensity of charged particles during a solar flare would essentially 
saturate the actual detector system and make the detection of prompt solar- 
cosmic-ray produced gamma-rays (or other gamma-rays) virtually impossible. 
Since the solar-cosmic-rays consist mainly of low energy protons with relatively 
few particles having energies above 100 MeV, the major induced activities in 
the moon are produced by (p,n), (p,pn), and similar low energy nuclear reac- 
tions. The dominant gamma-ray emitting isotopes produced by solar protons 
in the moon are Co 56 , Mn 54 , V 48 , Na 22 , and Al 26 . 

The production rates for solar -cosmic-ray induced radionuclides as a 

function of depth in the moon which were used in this work are those calculated 

by Reedy and Arnold (1972). These calculations were made assuming a rigidity 

spectral shape for the incident solar proton flux 

dJ 

— = kexp(-R/R 0 ) 


(4) 
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where R is the rigidity of the proton (R = pc/ze, p = momentum) in units of 
megavolts (MV) and R 0 is the parameter that determines the shape of the flux 
spectrum. Finkel et al (1971) found that production rates calculated using 
rigidity spectral shapes agreed very well with the observed activities of solar- 
cosmic-rays produced radionuclides in rock 12002, Details of the calculation 
of the production rates versus depth for solar-cosmic-ray induced radionuclides 
are given in Reedy and Arnold (1972). For all SCR-produced radionuclides, the 
production rates decrease rapidly with depth. 

Table 3 gives the calculated photon leakage intensities for the major solar- 
cosmic-ray induced radionuclides. Since the actual activity of each isotope de- 
pends on the intensity and spectral shape of the solar-cosmic-rays during the 
last few half-lives of the isotope (Reedy and Arnold, 1972), the results given 
here are for three different spectral shapes— values of R 0 of 50, 100, and 150 
MV— and for a steady state solar proton flux with the normalization that the 
integral omnidirectional flux of particles in space with energies above 10 MeV 
was 100 protons/cm 2 sec. These values are consistent with the solar-induced 
activities observed by Finkel et al (1971) and Wahlen et al (1972), which cor- 
responded to values of R 0 in the range of 80 to 100 MV and with equivalent 
steady state fluxes having integral normalizations above 10 MeV r angi ng from 
75 to 160 protons/cm 2 sec. The calculated intensities of Table 3 for Mn 54 , Na 22 , 
and Al 26 are about the same as the corresponding intensities resulting from the 
production of these nuclides from galactic -cosmic-rays (Table 4). Since their 
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half-lives are relatively long, yet shorter than the 11-year solar cycle, Mn 54 
and Na 22 have slowly varying solar-cosmic-ray induced activities during a solar 
cycle. The long-lived Al 26 has no variation in activity during a solar cycle. The 
very short-lived isotopes Co 56 (77-day) and V 48 (16-day) are made in only trace 
amounts by galactic -cosmic -rays and their activities will vary considerably 
during a solar cycle, depending on the solar flares during their last few mean 
lives. The agreement of the fluxes calculated here with those calculated by 
Armstrong (1972) is poor, due in part to the fact that he did not consider protons 
with energies of less than 30 MeV. 

For bodies further from the sun the SCR contribution will be even smaller, 
but the rate of decrease is not well known. 

Galactic-Cosmic-Rays . The galactic-cosmic -rays are high energy particles 
that continuously bombard the moon, producing many secondary particles as a 
result of the nuclear cascades they produce. The most important particle for 
inducing nuclear reactions in the moon is the neutron, approximately 10 neutrons 
being produced per incident galactic-cosmic-ray particle. These neutrons 
travel through the moon until they are removed by nuclear interactions or until 
they escape from the lunar surface. Secondary charged-particles are less im- 
portant since most of them have relatively low energies and usually are stopped 
by ionization-energy-losses before they have a chance to undergo nuclear 


interactions. 
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Neutrons with energies of about 1 MeV and greater can produce gamma- 
rays by exciting levels of nuclei by inelastic scattering reactions of the (n,x7) 


type, where x is usually a neutron but can be other particles. Galactic-cosmic - 
rays can also produce gamma-ray emitting radionuclides in spallation reactions. 
The production rates for both of these types of reactions were calculated using 
the galactic model of Reedy and Arnold (1972), Neutrons with energies below 
the first excited level of the nuclei in the moon can be elastically scattered by 
nuclei, escape from the moon, or be captured by nuclei through the (n, 7 ) 
reaction. In the process of neutron capture, prompt gamma-rays are emitted 
during the de-excitation of the product nucleus, which, in (n, 7 ) reactions is 
formed in excited levels that are about 8 MeV above their ground states. Some 
nuclei formed as a result of neutron capture are radioactive and emit gamma- 
rays in their decay. For the neutron capture reactions, the lunar neutron- 
capture rates of Lingenfelter et al (1961) and Lingenfelter et al (1972) were used,- 
as described below. 


The rates for inelastic scattering reactions and for the production of radi- 
onuclides at a depth x were calculated by integrating over all energies the 
product of the differential flux of galactic particles at that depth and the exci- 
tation function for the reaction 


OO 

A(x) = N f dE-^(E,x)o(E) 

Jo d E 


( 5 ) 
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with N being the number of target nuclei. The differential fluxes of galactic 
particles as a function of depth in the moon were taken from Reedy and Arnold 
(1972). For (n,X7) reactions, A(x) increases with depth until a depth of about 
50 g/cm 2 . Below that depth A(x) decreases with depth. For the production of 
radionuclides, the increase with depth is less, the reactions with the higher 
energy thresholds having less increase. 

Since the differential flux of galactic secondary particles decreases rapidly 
with energy (Reedy and Arnold, 1972) , most of the intense gamma-ray lines re- 
sult from reactions which have large cross-sections at low energies. Almost 
all of the (n,xy) reactions producing important gamma-ray lines have their 
largest cross-sections at relatively low energies, usually at energies below 
20 MeV and often as low as 1 MeV or less for some of the (n,ny) reactions. 

The reactions that produce the largest yields of radioactive nuclei are low- 
energy reactions such as the (n, p), (n,a), and (n,2n) reactions. In some cases, 
these radioactive nuclei decay through the same levels of their residual nuclei 
as are excited by (n,n7) inelastic scattering reactions on these nuclei. Thus 
the Al 27 (n,a)Na 24 and the Mg 24 (n, p)Na 24 reactions result in the same lines at 
1.369 and 2. 754 MeV as are produced by Mg 24 (n,n"y)Mg 24 reactions. Most of 
the excitation functions used in these calculations were based on experimental 


cross-sections. The references are cited in Table 4. 
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The calculated gamma- ray line intensities for the most intense lines for 
O, Mg, Al, Si, Ca, Ti, and Fe are given in Table 4 by element in order of de- 
creasing intensity. The intensities for gamma-ray lines resulting from neutron 
capture reactions are given and discussed separately below. Generally, the 
most intense line for each element for reactions other than neutron capture is 
that resulting from (n,n7) reactions exciting the first excited level of the most 
abundant isotope. Usually it is at an energy of less than 2 MeV. Several ele- 
ments have their first excited levels at higher energies, such asCa 40 and O 16 at 
3.73 and 6.129 MeV, respectively. Several elements have strong gamma-ray 
lines re stilting from prompt reactions other than the (n,n 7 ) reaction, such as 
the O l6 {n,a 7 ) C 13 lines between 3 and 4 MeV. 

The rates for the capture of low energy (thermal and epithermal) neutrons 
in the moon have been calculated for several chemical compositions by 
Lingenfelter et al (1961) and Lingenfelter et al (1972). Assuming a source 
spectrum of neutrons with energies below 10 MeV, the authors solved a multi- 
group transport equation to get the equilibrium energy spectrum and spatial 
distribution of the lunar neutron flux. They found that the calculated neutron 
flux was strongly dependent on the assumed composition of the lunar surface, 
the important compositional variables being the macroscopic capture cross- 
sections for the major elements, hydrogen content of the lunar surface, and the 


concentration of the rare earth elements Gd and Sm. 
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The presence of significant amounts of hydrogen in any region of the lunar 
surface would decrease the neutron leakage flux, increase the ratio of capture 
at low neutron energies to that at higher energies, and increase the capture rate 
in the top 150 g/cm 2 of the lunar surface (Lingenfelter et al, 1961). The total 
macroscopic capture cross-section (for nuclei with 1/v capture cross-sections) 
affects the energy spectrum since it changes the probability of being captured 
before reaching thermal energies if this cross section is larger. The rare 
earth el ;ments Gd and Sm have very large capture cross-sections near thermal 
energies and can capture a. significant fraction of the thermal neutrons. 

The total neutron capture rates of Lingenfelter et al (1972) and R. E. 
Lingenfelter and E. H. Canfield (private communication, 1971) were used in the 
calculation of the intensities of gamma-ray lines resulting from neutron capture. 
The rate for neutron capture increases significantly with depth to a depth of about 
150 g/cm 2 , below which the rate decreases with depth. The ratio of the neutrons 
captured to the neutrons escaping from the moon was not very sensitive to the 
total macroscopic neutron capture cross-section. The total macroscopic neutron 
capture cross-sections for the Apollo 12 lunar fines, the Apollo 14 lunar fines, 
and the lunar anorthosite are 0.78, 0.68 and 0.37, respectively, for that for 
the Apollo 11 lunar fines, but their total neutron capture rates were 0.97, 0.95, 
and 0.83, of the rate for the Apollo 11 fines. The main factor in determining the 
rate of neutron capture by a specific element is the fraction of the total macro- 
scopic neutron capture cross-section that the element contributes for that 
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chemical composition. The ratio of the fraction of the neutrons captured by an 
element to the weight fraction of that element is inversely proportional to the 
total macroscopic capture cross-section for a given chemical composition. 

For example, in the Apollo 11 lunar fines, iron (12.3% by weight) captures 35% 
of the neutrons, but in a lunar anorthosite (2.7% iron by weight), it captures 21% 
of the neutrons, whereas silicon (about 21% by weight in both) captures 8% and 
22% of the neutrons in these two chemical compositions, respectively. All 
elements except the rare earth elements were assumed to capture neutrons with 
1/v cross-sections. The contribution of the rare earth elements Gd and Sm were 
estimated using the equivalent 1/v cross-sections calculated by Lingenfelter 
et al (1972) for these elements. (For the Apollo 11 lunar fines, the effective 
cross- section for Gd is 18% of its thermal cross-section andGd captures 6% 
of all the neutrons.) The 1/v cross-sections used for the major elements were 
measured at 0.023 eV (2200 m/sec) and were taken from Bartholomew et al 
(1967), except for the value of 0.176 barns which was used for Si (Maddison and 
Jones, 1969), 

The yields per neutron capture for the prompt gamma-rays were taken from 
Bartholomew et al (1967). The neutron- capture gamma-ray energies and yields 
for neutrons captured above thermal energies were assumed to be the same as 
those for thermal neutrons, even though a major portion of the lunar neutrons 
are captured before they reach thermal energies. The yields per disintegration 
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for gamma-rays resulting from the decay of radioactive product nuclei were 
those of Lederer et al ,(1967) . The photon fluxes for the major lines of each 
element were calculated by numerically integrating Equation (2) and are given 
in Table 5 for the Apollo 11 chemical composition. 

The ratio of the gamma-ray flux for a given line to the weight fraction of 
the element varies with chemical composition. For the Apollo 12, Apollo 14, 
lunar anorthosite, and chondritic chemical compositions, this ratio of flux to 
weight fraction is 1.25, 1.40, 2.27, and 1.46 times that of the Apollo 11 compo- 
sition. 

The only difference among the chemical compositions affecting the neutron 
capture rates for all of these compositions but the chondritic composition was the 
total macroscopic cross-section and the fraction of the total captures each ele- 
ment contributed. The major difference between the Apollo 11 lunar fines and the 
chondritic composition was that the chondritic composition was arbitrarily as- 
sumed to include 0.1 atoms of hydrogen for each atom of silicon. This addition of 
hydrogen raised the ratio of the photon flux to the weight fraction for each element 
by a factor of about 1.4. A hydrogen to silicon atomic ratio of 1 would raise this 
ratio by a factor of about 2.7 above the ratio for no hydrogen. 

Gamma-ray lines for the rare earth elements Gd and Sm are not given here 
since they have few neutron capture lines with yields per capture greater than 
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a few percent. The main effect of these rare earth elements is to decrease the 
intensities of the major lines for all other elements. If there were no rare 
earth elements in the Apollo 11, 12, and 14 lunar fines, the intensities of the 
lines for all other elements would be raised by factors of 1.06, 1.11, and 1.25, 
respectively. Only for the Apollo 14 chemistry, where the elements with high 
thermal capture cross-sections, iron and titanium, are low and where the rare 
earth elements are very abundant, does the capture of neutrons by rare earth 
elements significantly affect other gamma-ray line intensities. 

The calculated production rates for (n,xy) and (n,y) reactions are for the 
average GCR flux near the earth during a solar cycle. Since the GCR flux is 
somewhat lower at solar maximum and higher at solar minimum than the aver- 
age flux used here (Webber, 1967) , the actual production rates will vary with 
the solar cycle. Thus the gamma-ray flux for all sources except those from 
the decay of long-lived isotopes will vary during the solar cycle. Since all GCR 
produced lines will have approximately the same variations, only the average 
flux for the gamma-ray lines are given here. A modest scale correction de- 
termined from cosmic -ray data must be applied for the time of each mission. 
The gamma-ray fluxes expected from other large bodies in the solar system 
should be similar to those calculated here, as the intensity of the high energy 
GCR particles does not vary significantly with distance from the sun. 
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The fluxes calculated for the prompt GCR produced lines agree fairly well 
with those calculated by Armstrong (1972) using a Monte Carlo method. The 
currents for the neutron capture lines agree well with the calculated currents of 
Gorenstein and Gursky (1970), while their results for the lines resulting from 
(n,x7) reactions are about twice our values. The disagreement for the (n,xy) 
reactions is probably due to the fact that they used a simpler model for these 
reactions than the one used here. 

Calculated Fluxes. The fluxes from a semi-infinite plane of the major 
gamma-ray lines for each of the five chemical compositions of Table 2 are 
shown in Figures 1 and 2. The fluxes include the contributions from the decay 
of the natural radioactive nuclides, from all the galactic-cosmic-ray sources, 
including both prompt-(n, x7)and (n, 7 ) — and decay photons, and from the decay 
of solar-cosmic-ray produced nuclides for an incident solar proton spectral 
shape R 0 = 100 MV and a 4 it proton integral flux above 10 MeV of 100 protons/ 
cm 2 sec. The line at 511 keV resulting from the annihilation of positrons was 
not calculated, because of the many different sources for this line. It is simply 
indicated in the figures. In addition to the lines shown in Figures 1 and 2, there 
will be a strong continuum that includes a multitude of minor lines, photons re- 
sulting from Compton scattering and bremsstrahlung in the moon. The approxi- 
mate nature of this continuum is discussed below. 

The major lines shown in Figures 1 and 2 are also given in Tables 1, 3, 4, 
and 5 for the Apollo 11 chemical composition. There are only a few gamma-ray 
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lines that stand out above the rest. Most of the major lines have about the same 
magnitudes and there are a large number of such lines, especially below 4 MeV. 
The major elements whose abundances in the moon are expected to be similar 
for all likely chemical compositions— oxygen and silicon— have many intense 
lines. Several elements, such as magnesium, aluminum, and calcium, have 
most of their most intense lines near other major lines. For example, the 
Mg 24 (n, ny) line at 1.369 MeV is at an energy where there are many other lines 
from such sources as U, Na 24 , and Ti(n, 7 ). The Ca(n,n 7 ) lines at 3.73 and 
3. 90 MeV are near the 0 16 (n,a7) lines at 3.68 and 3. 85 MeV. 

The calculations indicate that there are good gamma-ray lines with which 
to determine the concentrations in the moon's surface of oxygen, iron, and 
titanium. Concentrations of thorium, potassium, and uranium can also be 
readily observed. Calcium, silicon, magnesium and aluminum also have some 
good gamma-ray lines, although few of their lines are well isolated from other 
major lines, and, while absolute concentration will be hard to determine, large 
variations of these elements in the moon should be observable. The other ele- 
ments do not have very useful gamma-ray lines for expected range of 
compositions . 

SPATIAL RESOLUTION 

To understand the possible spatial resolution of the lunar experiment, it is 
necessary to consider the angular dependence of the system response. Figure 3 
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shows the basic geometry. For a given detector-source distance £ the photons 
striking the detector are emitted at a fixed angle 0 from the normal at the source, 
and strike the detector at an angle 6 with respect to the line connecting the de- 
tector and the center of the moon. The limiting distance, £ max , is the horizon 
distance, equal to (2Rh + h 2 ) ,/2 . For a height of 100 km above the moon's sur- 
face, £ max is 598 km. 



By integrating over radius the product of Equation (6) and the moon's sur- 
face area at each radius, one can calculate the flux at the detector originating 
within a circle of radius p on the moon's surface. The values for a , J3 , and £ 
in Equation (6) all vary withp . It is of interest also to consider the behavior 
in crossing a great circle boundary between two extended source regions of 
different composition, at right angles to the boundary. The flux for beyond the 
boundary is calculated the same as for the flux for inside a circle, with the 
integration over radius including an additional factor for the fraction of the 
circle at radius p beyond the boundary. 
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These calculations were done for four cases: an (n, 7 ) reaction (the 4.936 
MeV line from Si(n, 7 )), an (n.xy) reaction (the 1.369 MeV line from Mgfn.ny)), 
a source uniform with depth (any natural radioactive decay line), and a solar- 
cosmic-ray produced radionuclide (the 1.809 MeV line of Al 26 produced by 
protons on Si). The difference among these cases is the behavior of A(x) in the 

first few photon mean free paths of the surface: A(x) increases considerably with 

the depth near the surface for the (n,y) reaction. Increases seme what with depth 

for the (n, x7 ) reaction, and decreases rapidly with depth for the solar-cosmic- 
ray (SCR) case. The relative fractions of the flux at an isotropic detector 100 km 
above the lunar surface are given in Figure 3 for point sources, inside a circle, 
and beyond a boundary, for each of these four cases. 

For the uniform source and a detector altitude of 100 km, the distance p 
where F falls to half the value of the point source for directly below the detector 
is 73 km. Half the flux at the detector originates inside a circle of radius 118 
km. At 68 km from a great circle boundary 75% of the field of view is occupied 
by the region below the detector and 25% by the region beyond the boundary. 
Corresponding distances for the SCR case are considerably greater than these. 

For the (n,x 7 ) and (n, 7 ) cases, these distances are about 90% and 80%, re- 
spectively, of those of the uniform case. GCR-produced radionuclides vary be- 
tween the uniform case and the (n,x7) case, depending on the excitation function 


for the reaction. 
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The SCR-produced radionuclides, where A(x) decreases with depth, are 
examples of limb brightening, where more photons arrive per unit solid angle 
at the detector from regions near the limb. The (n,7) and <n,x7) cases, where 
A(x) increases with depth near the surface, are examples of limb darkening. 
Since the flux of photons per unit solid angle at the detector varies with angle 
for non-uniform sources, the flux at the detector above the moon is not just 
the flux for a semi -infinite plane times the ratio of the solid angle for the lunar 
case to that for the semi -infinite plane. For an altitude of 100 km above the 
moon, this ratio of solid angles is 0. 67465, 0 max being 71.013°. For this lunar 
altitude, the ratio of fluxes from the moon to that from a semi-infinite plane is 
greater than the corresponding ratio for solid angles by a factor of about 1.12 
for most (n,7) cases, 1.07 for most (n,x7) cases, and less by factors of about 
0.7 to 0.8 for various SCR eases. For sources uniform with depth, the ratio of 
fluxes is exactly the same as the ratio of solid angles for all geometries. 


DATA REDUCTION 

As discussed above the lunar gamma-ray emission spectrum can be considered 
to be made up of two components: the line spectrum and a spectrum con- 
tinuous in energy. This continuum can be associated with such processes as 
Compton scattering of gamma-rays, and bremsstrahlung in the lunar surface. 

There are also a number of other contributions to the background continuum, 
such as induced and natural radioactivity in the spacecraft, diffuse sky back- 
ground, and spallation effects in the detector crystal. From laboratory 
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experiences and from the analysis of earlier space flight gamma-ray data 
we find that the total background continuum can be approximated quite closely by 
a power law distribution in energy. In this section we will describe the analytic 
method for obtaining the elemental composition from the observed lunar surface 
spectrum. 

We first consider the basic method of inferring photon spectra from an 
analysis of the measured pulse height spectrum. The pulse height spectrum is 
a measure of the interaction between the incident gamma-ray and the detector. 
We can describe this interaction process as follows: 

rE 

YfV) = I maX T(E) S(E,V)dE (8) 

J 0 

where T(E) is the differential energy photon spectrum as a function of energy E, 
S(E, V) is the detector response function which reflects the mechanisms of con- 
verting photons of energy, E, into a pulse height signal, V, at the output of the 
detector, and Y(V) is the measured "pulse height" spectrum as a function of 
pulse height V. 

The problem is to infer T(E) from Y(V) and a knowledge of the function 
S(E, V). We now consider the shape of S(E, V) and its variation with energy. 

S(E, V) can be either experimentally determined or theoretically calculated, 
(Bell, 1955; Seltzer and Berger, 1971). Radioactive sources were available 
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to study this energy dependence in the Apollo detector system (a 2-3/4 x 2-3/4 
inch Nal(TC) crystal) experimentally up to 3 MeV. Neutron capture spectra 
were used to study crystal response up to 10 MeV (Trombka et al, 1971), The 
energy response functions were also calculated using Monte Carlo techniques and 
found to be in good agreement with experiment (Trombka et al, 1971; Seltzer and 
Berger, 1971). We conclude that the response function can be calculated 
theoretically for any energy needed in the analysis considered in the following 
discussirn. The three pulse height spectra shown in Figure 4 have been produced 
by the Monte Carlo calculation technique. 

The analysis of gamma-ray pulse height spectra in this spectral region must 
be accomplished by numerical methods because we cannot yet perform the in- 
verse transform analytically. The method for obtaining the best estimate for 
T(E) is based on a theorem in sampling theory; and on the fact that there is a 
finite energy resolution of the detection system. The Shannon Sampling Theorem 
asserts that, if a distribution has no oscillatory component with a frequency greater 
than f max , then samples at discrete points not farther apart than 1/2 f max 
describe the original function exactly. The original function may then be re- 
constructed from the samples (Trombka et al, 1970). The approximation formula 
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replaces the integral in Equation (8) with a sum j 

Y. = 2 (9) 

j J J 

where Yj is the value of Y(V) corresponding to channel or pulse height i, Tj is 
the total number of gamma-rays in the energy interval AEj about the energy Ej , 
and Sy is the intensity in channel i for a gamma-ray of the jth energy group. In 
applying Shannon's Sampling Theorem to the analysis of spectra obtained with 
Nal(Tfi) crystals, sampling components are chosen so that the energy intervals 
between the energy groups are less than the half width at half maximum of the 
photopeak at that energy (Linden, 1959; Trombka and Schmadebeck, 1968). 

A valid numerical transformation can now be obtained using techniques de- 
rived from the analysis of variance. The details of this derivation can be found 
in Trombka et al (1970). Using the least square principle, the method requires 
that 

M - 2 ojj ( Yj - 2 TjSj:) 2 (10) 

i j 

be a minimum, where the statistical weight corresponding to the measurement 
of Yj in channel i, ojj = 1 /of (where of is the variance in the measurement of 
the count in channel i). The variance of data processing must include all effects 
such as background subtraction. The other terms have been defined. The 
partial derivatives of M are taken with respect to T in order to determine the 
minimum. The solution in matrix form is given by 


T = (SwS)' 1 SojY 


(ID 
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where Y is a vector describing the pulse height spectrum, S is an m x n matrix 
describing the discrete set of detector response function, S is the transposition 
of S, oj is a square diagonal matrix of the weighting functions, and T is a best 
estimate of the differential energy photon spectrum. A detailed description of 
the calculating method and computer program is found in Trombka and 
Schmadebeck (1968). 

We are now able to transform from measurement space to photon space. A 
major problem exists because we have a mixture of discrete lines and continuum. 
This problem is further complicated in the determinations of the inverse matrix 
(ScoS) -1 . One obtains oscillations in the solution due to interference between 
the components of the S matrix. These oscillations become most pronounced in 
the vicinity of discrete lines. Methods for damping out these oscillations have 
been developed. Our method essentially consists in trying to obtain a descrip- 
tion of the continuous component using the transposition from pulse height to 
energy space (Equation 11), eliminating all possible discrete lines and points on 
either side of those lines where oscillation occur, and finding the best fit of a 
simple spectrum, generally a power law distribution, to the remaining points. 

Our analytic methods allows the use of other types of analytic and empirical 
background shapes. In fact if one wishes, a number of interations on the analysis 
can be performed using various background shapes in order to find the form 
most suitable for the particular situation. For purposes of illustration of the 
analytical method we will assume the background can be described by a power 



law distribution. We thus assume that the continuum can be approximated by 
the function 

T c = AE -n (12) 

where A is an amplitude constant, and n is the power to be determined by fitting 
those points which characterize the continuous distribution. 

Once T c is determined, the corresponding pulse height distribution can be 
obtained by using Equation (9). The pulse height spectrum of T c is then sub- 
tracted from the measured spectrum, and the pulse height distribution of the dis- 
crete component obtained. This can then be converted to a photon spectrum. The 
geochemical information can then be deduced from this final photon distribution. 

An example of the use of the analytical method may be a help in understanding 
its application. Using the methods of the first part of this paper the discrete line 
gamma- ray emission spectrum from the lunar surface for a given composition can 
be obtained. An average lunar surface composition for the major elements com- 
position was assumed and is indicated in Table 6. The expected photon flux is 
also given. For the neutron capture lines, the flux to weight ratios used were 
1.4 of the Apollo 11 values. Using Monte Carlo methods, a pulse height spectrum 
due to the discrete lines observed at the orbital altitude with Nal(TC) detector 
was generated. The detector resolution for 0.661 MeV was a 2-3/4" by 2-3/4" 
assumed to be 8. 5%. We now assumed that the continuum spectrum was of the 
form T c = AE' L3 . From the Luna 10 data (Vinogradov et al, 1967) we estimated 
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that the continuous distribution was about 85% of the total observed flux from 
the lunar surface. Finally, it was assumed that the spectrum was accumulated 
for 200 minutes. The calculated spectrum is the top curve shown in Figure 5, 

The spectrum was divided into 511 channels. Random fluctuations due to counting 
statistics were included. Instrument parameters correspond to the system used 
during the Apollo mission. 

This pulse height spectrum was transformed back to a photon spectrum 
using a litrary of 118 monoenergetic components. The results are displayed 
in Figure 6. The oscillations around lines and at high energies can be easily 
seen. These oscillations and discrete lines were eliminated. Seventeen points 
remain to characterize the continuous distribution (the points circled). A least 
square fit was performed assuming a power law dependence. Both the magnitude 
A and power n were determined. For this example it was found that A = 0.062 
and that n = 1. 33, only 2. 3% from the value assumed in the calculation. The 
calculated power law distribution of the continuum is shown as the line in Figure 
6 and is assumed to extend up to 10 MeV. Using the values of A and n obtained 
from the least square fit, the shape and magnitude of the pulse height spectrum 
of the continuous spectrum can be obtained. It is the smooth middle curve of 
Figure 5. This component is subtracted from the top spectrum shown in Figure 
5, and an estimate of the discrete line spectrum is obtained. This is the bottom 


curve of Figure 5. 
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To analyze this line spectrum, pulse height spectra characteristic of the 
elemental reaction components (rather than monoenergetic gamma-ray lines) are 
used for the library. All gamma-ray lines for natural radioactive decay, neutron 
capture, and neutron inelastic scatter for each element are grouped separately. 
For example the pulse height spectrum for the natural emission spectrum from 
Th and that expected from the O (n,x 7 ) reactions are shown in Figure 7. The 
set of components used to form a fifteen component library is shown in Table 
6. Using this library, the elemental photon flux and elemental composition can 
be obtained by matrix inversion (Equation 11). 

The results of the analysis of the discrete pulse height spectrum are given 
in Table 6. These recovered results can he compared with the values assumed 
in the calculation. If our assumptions about the continuum and discrete line 
spectra are correct, our capabilities for obtaining qualitative and quantitative 
elemental analysis from the measured gamma-ray are defined by these results 
and are discussed below. 

The errors indicated in the table are associated with counting statistics, 
background subtraction, and mutual interference of components. Our ability to 
determine the natural radioactive elements, Th, U, and K should be limited 
only by counting statistics. The results agree well with the true values. The 
reactions Mg(n,n 7 ), Fe(n,n7), and Si(n,ny) can he used to determine the 


elemental weight fraction of Mg, Fe, and Si. 
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One of the major sources of error in the determination of the elemental 
weight fraction can be attributed to interference between various spectral shapes, 
that is the inability to resolve elemental components separately. An analysis of 
the co-variance between various elemental components making up the pulse 
height spectra of the libraries used in the analysis of the discrete spectra has 
been carried out and the percent interference between various components cal- 
culated. The method of calculation is described in Trombka and Schmadebeck 
(1968). The higher the percent interference, the more difficult is it to resolve 
the particular components. Our laboratory experience indicates that inter- 
ferences greater than 20% cause trouble in interpretation. Table 7 indicates 
the percent interferences greater than 20% between the various library com- 
ponents used in the analysis of the elemental reaction pulse height spectra. For 
example one finds that the Fe(n, 7 ) spectrum strongly interferes with the 
0(n,x7), and Si(n,n 7 ) and Al(n, 7 ) spectra. Thus it is hard to use these re- 
sults to determine Fe concentration. This analysis also predicts difficulties in 
using other high energy lines. Some of the problem can be removed if it is pos- 
sible to use plausible geochemical assumptions. For example, it might be 
reasonable to assign the concentration of Fe from the Fe(n,n 7 ) component, or 
of Si and O from general geochemical considerations. 
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Chi-square is also calculated to determine the goodness of fit. This param- 
eter is of utmost interest in determining if there has been any significant change 
in the gain or zero level of the spectrometer, and whether a component or com- 
ponents are missing. For details on the application of these parameters see Trombka 
and Schmadebeck ( 1968 ) and Trombka et al (1970 ). Great care must be used 

in interpreting the Chi-square fit. A good fit does not necessary represent truth 
in terms of geochemical analysis. Thus using monoelemental library components 
in the least square fit, yield poorer Chi-square values than using monoenergetic 
library components but one obtains more physically meaningful results with 
the elemental library components. 

The procedure described above illustrates one iteration on the data. In 
actual practice a number of iterations are used. For the second iteration one 
eliminates all monoelemental components with negative intensities or with 
intensities with an error greater than 100%. The remaining elemental components, 
which correspond for the most with the higher concentration elements, are then 
used to construct a composite discrete energy pulse height spectrum. This 
discrete pulse height spectrum is then subtracted from the original and a better 
estimate of the background continuum is obtained. The analysis as described 
above can then be repeated. 

The procedure for analyzing the measured pulse height spectra described 
above allows us to perform analysis of the Apollo 15 and 16 gamma-ray spectra. 
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Both qualitative and quantitative results can be obtained. Accuracy and pre- 
cision of results are determined in terms of counting statistics, background 
subtraction, and spectral interference. 


DISCUSSION 

Although we have discussed the use of ’'ground truth" above, we must re- 
turn to it at this stage. The matrix inversion program can make some use of 
ground truth; for example, the ratio Th/U may be set at a geochemically known 
value around 3,8 or the concentration of O and Si may be assumed constant in 
starting an iteration. Beyond this, the final chemical concentrations derived 
may be compared with the soil analyses in suitable regions such as Mare 
Tranquillitatis. This has been and will be useful in developing the computer pro- 
cedures to higher precision and accuracy. It may also be necessary at a future 
stage to introduce a scale factor for one or more elements whose analysis de- 
pends on cosmic-ray interaction, either because of deficiencies in the Reedy- 
Arnold model or limitations of our knowledge of cross-sections. Other sorts of 
experimental evidence may enter at this point. Simulation experiments using 
accelerators, and balloon flights which carry a detect© r-target package to the 
top of the earth's atmosphere are both important here. 

While the general sky background seems to consist mainly of continuum, 
and hence is included in the empirically derived lunar continuum, its study is 
necessary as a check of our understanding. The spacecraft contributes a small 
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but measurable line spectrum. The lines of K and Th come from known sources 
and have been measured on the ground. Other lines are produced by cosmic- 
ray bombardment. All these are best extracted from the spectrum observed in 
transearth coast, since the continuum there is much lower in intensity than the 
lunar one. The corrections for K and Th are significant at the lowest levels 
observed on the moon. 

At present the list of elements which can be determined with reasonable 
accuracy on the moon is as follows: Th, U, K, Si, O, Fe, Mg. As methods 
are refined the data already obtained are expected to yield more (Metzger, 

1972). 

For future missions, the use of a Ge(Li) detector system, possibly with 
collimation to improve areal resolution, would be a major improvement. The 
much higher energy resolution would greatly simplify the required data analysis. 

In extending the process to other objects in the solar system, there are a 
few general guidelines. First, while the minimum necessary counting time for 
meaningful results on major elements will depend on closeness of approach, 
size of object, instrument design, and composition, the value observed so far 
in our Apollo analyses, of the order of 30 minutes, is encouraging for flyby 
missions to larger objects. The minimum time required for detailed mapping 
of an entire object on a resolution of the order of 5° will be of the order of a 
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month. Since the K/U ratio on the lunar surface is lower than that observed in 
meteorites or on earth, we might anticipate that the spectrum of another rocky 
object would show the K 40 line at 1.46 MeV more prominently. The spectrum of 
an object much higher in hydrogen, such as a comet, will be qualitatively dif- 
ferent, with neutron capture lines much more prominent. Finally, Mars is a 
special case, with an atmosphere comparable to the mean thickness of the lower- 
energy gamma-ray lines. The experiment is still quite feasible; the special 
aspects are discussed elsewhere (Metzger et al, 1970; Surkov et al, 1971) 
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Table 1 

Gamma-ray fluxes from the decay of natural radioactivities. The half- 
lives used for K 40 , Th 232 , and U 238 were 1. 26 x 10 9 , 1.41 x 10 10 , and 
4.51 x 10 9 years, respectively, and the K 40 /K ratio used was 0.000118. The 
fluxes are for the Apollo 11 chemical composition (K, Th, and U values 


of 1150, 2. 1, and 0. 55 ppm by weight, respectively). 


Parent 

Element 

Isotope 

Energy 

(MeV) 

Yield/ 

disintegration 

Flux (photons/ 
cm 2 min. ) 

K 

K 40 

1.46 

0.11 

2.37 

Th 

oo 

O 

rN 

i— i 

2.614 

0.36 

2.41 


Ac 228 

1.593 

0.045 

0.24 


Ac 228 

0.97 

0. 20 

0. 82 


Ac 228 

0.91 

0. 23 

0. 92 


Bi 212 

0. 727 

0.07 

0. 25 


rp|208 

0.583 

0.31 

1.02 


ip^208 

0.511 

0. 083 

0. 26 


Ac 228 

0.337 

0. 15 

0.39 


Pb 212 

0, 239 

0.47 

1.06 

U 

Bi 214 

2,204 

0.05 

0.24 


Bi 214 

1.764 

0.17 

0.74 


Bi 214 

1.378 

0.05 

0. 19 


Bi 214 

1.238 

0. 06 

0.22 
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Table 1 (Continued) 


Parent 

Element 

Isotope 

Energy 

(MeV) 

Yield/ 

disintegration 

Flux (photons/ 
cm 2 min. ) 


Bi 214 

1.12 

0.17 

0. 59 


Bi 214 

0. 769 

0.05 

0. 15 


Bi 214 

0.609 

0.47 

1. 24 


Pb 214 

0.352 

0.36 

0. 76 


Pb 214 


0.295 


0.19 


0.37 
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Table 2 


Chemical compositions used in calculations. Abundances of 
elements are in percent or parts per million (ppm) by weight. 



Apollo 

11 

Apollo 

12 

Apollo 

14 

Lunar 

Anorthosite 

Chondrite 

o 

41.6 

42.6 

44.2 

45.0 

35.0 

Na 

0.32 

0.32 

0.47 

0.3 

0.68 

Mg 

4.6 

6.3 

5.6 

2.4 

14.4 

A1 

7.1 

6.8 

9.2 

16.0 

1.3 

Si 

20.0 

21.4 

22.5 

20.7 

18.0 

K 

0. 115 

0.18 

0.43 

~0.04 

0.086 

Ca 

8.6 

7.2 

7.6 

11.5 

1.4 

Ti 

4.6 

1. 7 

1.0 

0.2 

0.08 

Fe 

12.3 

12.0 

8.0 

2.7 

25.0 

P(ppm) 

600 

1200 

2100 

400 

1000 

S(ppm) 

1000 

— 

— 

— 

23000 

Cr(ppm) 

1950 

2900 

1300 

— 

2500 

Mn{ppm) 

1600 

1700 

1000 

— 

2600 

Ni(ppm) 

240 

160 

— 

— 

13000 

Gd(ppm) 

17 

20 

35 

~1 

0.34 

Th(ppm) 

2.1 

8 

14 

-0.7 

0.04 

U(ppm) 

0. 55 

2 

4 

~0.2 

0.013 

Sum 

99. 78 

99. 10 

99.45 

98.88 

100. 16 
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Table 3 

Gamma-ray fluxes from solar -cosmic -ray produced radionuclides. The 
results are given for three different incident solar proton exponential 
rigidity spectral shapes and for 4^ integral proton fluxes above 10 MeV 
of 100 protons/cm 2 sec. The gamma-ray fluxes are for the 
Apollo 11 chemical composition. 

Fluxes {Photons/ cm min) 

Target Radionuclide Energy 

(MeV) R 0 = 50 MV 100 MV 150 MV 


Na 

Na 22 

1.275 

0.008 

0. 022 

0.032 

Mg 

Na 22 

1.275 

0.024 

0. 102 

0.189 

A1 

Na 22 

1.275 

0.016 

0.084 

0. 152 


Al 26 

1. 809 

0. 149 

0.360 

0. 518 

Si 

Na 22 

1.275 

0.012 

0.104 

0. 241 


Al 26 

1. 809 

0.078 

0.378 

0.689 

Ti 

y48 

0.983 

0.183 

0.175 

0.162 


V 48 

1.312 

0.183 

0.177 

0.165 

Fe 

Mn 54 

0. 835 

0.033 

0.140 

0. 240 


Co 56 

0. 847 

0. 340 

0. 334 

0.306 


Co 56 

1,238 

0.235 

0. 233 

0. 215 


50 


Table 4 

Gamma-ray fluxes for lines resulting from galactic-cosmic-ray 
induced inelastic scattering and radioactivities. References are the 
sources for the cross-sections used in calculating the production 
rates for these reactions. The fluxes are for the Apollo 11 
chemical composition and are in order of decreasing intensity 


for each element. 


Reaction 

Energy (MeV) 

Flux (Photons/cm 2 min) 

O l 6 (n,nT) (a) 

6. 129 

2.43 

0 16 (n, n 7 ) (a) 

7. 117 

0. 94 

O 16 (n, n7) (a) 

6.917 

0. 79 

O 16 (n,a 7 ) (a) 

3.68 

0.56 

0 1 <s <n,n7) (a) 

2.743 

0.46 

0 16 (n,na 7 ) (a) 

4. 433 

0.44 

0 16 (n,a 7 ) (a) 

3. 85 

0.28 

O 16 (n,p)N 16(b) 

6.129 

0 . 26 

0 16 (n,a 7 ) (a) 

3.09 

0.19 

Mg 24 (n, n 7 ) (b) 

1.369 

0. 78 

Mg 24 (n,p)Na 24(b> 

2.754 

0. 12 

Al 27 (n,n 7 ) (c) 

2.21 

0.43 

Al 27 (n,n 7 ) (c) 

1. 013 

0. 40 

Al 27 {n,n7) (c) 

3.001 

0.28 



51 


Table 4 (Continued) 


Reaction 

Energy (MeV) 

Flux (Photons/cm 2 min) 

Al 27 (n, 2n)Al 26<d) 

1.809 

0. 19 

Al 27 (n, n-y) {c) 

0. 843 

0. 16 

Al 27 (n,a )Na 24(b) 

2. 754 

0. 13 

Si 28 (n, n7) (e) 

1. 779 

2.65 

Si 28 (n,p)Al 28(b) 

1. 779 

0. 72 

Si 28 (n,pa)Na 24(0 

2. 754 

0.31 

Si 28 (n, p2n)Al 26(d) 

1. 809 

0.31 

Si 28 (n, n 7) (e) 

2.835 

0.27 

Ca 40 (n,n 7 ) (8> 

3. 73 

0.22 

Ca 40 (n,n 7 ) (s) 

3. 90 

0.13 

Ti 48 (n,n7) (f) 

0. 983 

0.29 

Fe 56 (n,nT) <c) 

0.847 

1.53 

Fe 56 (n,n7) (c> 

1. 238 

0.37 

Fe 56 (n, n 7 ) (c) 

1.811 

0.16 

Fe 56 (n, n 7 ) (c) 

2. 113 

0.12 

Fe(n, x)Mn 54(d> 

0.835 

0.11 


^Orphan, Hoot, and John (1970) 
(b) Stehn et al (1964) 

(^Orphan and Hoot (1971) 
^^Reedy and Arnold (1972) 



Table 4 (Continued) 

( e )Dickens (1970) for neutron energies below 10 MeV and Drake (1968) for 
neutron energies from 10 to 20 MeV 
(0 Cros s -sections estimated from analogous reactions 
^Drake (1967) 
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Ca(n, 7 ) 


Ti(n,7 ) 


Fe(n,7) 


.92 
7.724 
1.779 
4. 936 
3.54 
6.42 
1.943 
6. 755 
6.42 
1.38 
9.298 
7.643 
7.629 
7.277 
6.018 


0.36 
0.23 
1.00 
0.53 
0.53 
0.40 
0. 89 
0.47 
0.31 
0.86 
0.03 
0.21 
0. 21 
0.05 
0.08 
0.08 


0.02 
0.09 
0.14 
0.34 
0.27 
0.22 
0.21 
1.67 
1.07 
0.97 
0. 12 
0. 75 
0. 75 
0.18 
0. 25 
0. 25 


5. 92 



Synthetic and 



0(n,xY) 

6.60 

Mg(n, n7) 

1.32 

Al(n, n 7 ) 

3.06 

Al(n, 7 ) 

0.70 

Fe(n, n 7 ) 

2. 23 

Fe(n, 7 ) 

3.10 

Si(n,n 7 ) 

5.49 

Si(n,7) 

1 . 27 

Ca(n,n 7 ) 

0.48 

Ca(n, 7 ) 

0.84 

Ti(n, n7) 

0. 16 


Ti(n, 7 ) 


1. 88 
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e 6 

id Weight Fractions and 

e Lunar Composition 

Recovered 

Weight Fraction Weight Fraction 


2. l(ppm) 

2. 11 

± 0.14 

0. 55 (ppm) 

0. 53 

± 0.04 

1150(ppm) 

1168 

± 92 

0.42 

0.49 

± 0.05 

0. 05 

0. 060 ± 0. 012 

0.10 

0. 15 

± 0.02 

0. 10 

0.15 

± 0.17 

0. 10 

0. 10 

± 0.01 

0. 10 

0.14 

± 0. 02 

0.21 

0.21 

± 0. 03 

0. 21 

0. 35 

± 0. 08 

0.09 

0. 22 

± 0.07 

0.09 

0. 11 

± 0. 05 

0.015 

0. 010 ± 0. 003 


0. 015 


0. Oil ± 0. 005 



Table 7 


Percent Interference Between Library Components 



Mg 

Th 

(n, n 7 ) 

U 

O 

(n,xy) 

A1 

(n.n7) 

Fe Si Ca 

<n,n7) <n,n7) (n, n7) 

Ti 

(n,n7) 

A1 

(n,T) 

Si Ca 

(n,7) (n,7) 

Ti Fe 

K 40 

(n» 7) (n, 7 ) 

Th 

100 

20.2 









Mg{n, nT) 

100 


21.0 




29.9 




U 

20.2 

100 









0(n,xy) 

Al(n,nT) 
Fe(n, ny) 

21.0 


100 

100 

100 

31.5 




51.2 22.4 

Si(n,n7) 






100 


73.6 


56.7 

Ca{n, ny) 
Ti(n,ny) 



31.5 



100 

100 



21.5 

Al(n, 7 ) 






73.6 


100 


79.2 

Si(n, 7 ) 









100 


Ca(n, 7 ) 









100 


Ti(n,7 ) 

29.9 


51.2 







100 

Fe(n, 7 ) 



22.4 



56.7 21.5 


79.2 


100 

K 40 










100 


cn 

cn 



FIGURE CAPTIONS 


Figure 1. The calculated gamma-ray fluxes for the indicated chemical compo- 
sitions at an isotropic detector just above the moon's surface (2 ir 
geometry). Except for the line at 0.511 MeV, the fluxes are for the 
discrete lines produced by the decay of naturally occurring radio- 
nuclides (Table 1) and from the interaction of GCR and SCR particles 
with the moon (Tables 3-5). The data are grouped in energy intervals 
50 keV wide. Note that the flux scale is the same as that in Fig. 2. 

Figure 2. The calculated gamma-ray fluxes for the indicated chemical compo- 
sitions at an isotropic detector just above the moon's surface (2 it 
geometry). Except for the line at 0.511 MeV, the fluxes are for the 
discrete lines produced by the decay of naturally occurring radio- 
nuclides (Table 1) and from the interaction of GCR and SCR particles 
with the moon (Tables 3-5). The data are grouped in energy inter- 
vals 50 keV wide. Note that the flux scale is the same as that in 
Fig. 1. 

Figure 3. The top diagram shows the angles (0 and 8) and distances for a 

detector at an altitude h above the moon relative to a point source on 
its surface. The graph shows the relative fluxes at an isotropic 
detector from a point source (I' # Equation 6), from inside a circle, 



and from beyond a great circle boundary, as a function of the 
distance along the surface of the moon from the sub-detector point. 

The fluxes are for an altitude of 100 km above the moon (B moon = 

© 

1738 km) and for four different profiles for gamma-ray source 
intensity versus depth. 

Figure 4. The pulse height spectra for monoenergetic gamma-rays at 3 dif- 
ferent energies incident upon a 3 by 3 inch NaI(T £) detector (~19 
keV/channel). The important spectral features are identified. 

These spectra are theoretical Monte Carlo calculations. 

Figure 5. The top curve is a pulse height spectrum calculated for the average 
lunar chemical composition given in Table 6 and includes both dis- 
crete lines and continuum as described in the text. This spectrum 
is for a 2-3/4 by 2-3/4 inch Nal(Tfi) crystal (~19 keV/channel) and 
includes statistical errors appropriate for a 200 min. integration 
time. The middle curve is the pulse height spectrum of the con- 
tinuum shown in Fig. 6. The bottom curve is the pulse height 
spectrum for the discrete lines obtained when the continuum con- 
tribution is subtracted from the original spectrum (top curve) . 

Figure 6. The dots are the lunar photon spectrum obtained when the synthesized 
lunar pulse height spectrum (top curve. Fig. 5) was transformed 




by a matrix inversion (Equation 11) using 118 monoenergetic com- 
ponents. The 17 circles are the intensities at energies where there 
are no discrete lines used to obtain the continuum and the straight 
line is the least squares fit to these points, the best fit being T c = 
0.062 E' 1,33 . 

Figure 7. The pulse height spectra used for the monoelemental components Th 
and 0(n,X7), calculated for a 2-3/4 by 2-3/4 inch NaI(T C) detector 
(~19 keV/channel). These components, along with those for the 
other elemental components given in Table 6, were used for the 
matrix inversion of the discrete pulse height spectrum (bottom 
curve, Fig. 5) in obtaining the results given in Table 6. 
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